Abnormal dilepton yield from local parity breaking in heavy-ion collisions 
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We propose a novel explanation for the dilepton excess observed in heavy ion collisions at invariant 
masses below 1 GeV. We argue that the presence of local parity breaking due to a time-dependent 
isosinglet and/or isotriplet pseudoscalar condensate substantially modifies the dispersion relation of 
photons and vector mesons propagating in such a medium, changing the p spectral function and 
resulting in a potentially large excess of e + e~ and p + pT with respect to the theoretical predictions 
based in a 'cocktail' of known hadronic processes. Possible signatures to prove or disprove this effect 
are discussed. 
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During the last decade several experiments in heavy 
ion collisions have indicated an abnormal yield of lepton 
pairs with invariant mass M < 1 GeV in the region of 
small rapidities and moderate transversal momenta [l], 0] 
(reviewed in [3J, |jj ) ■ This effect is visible only for colli- 
sions that are central or semi-central, and is present for 
e + e _ and also for p + p~ . From a comparison to pp and 
p-nuclcus collisions it has been established beyond doubt 
that such an enhancement is a nuclear medium effect Q . 
For the energies accessible at GSI (of few GeV per nu- 
cleon in HADES experiments |3[) the effect was partially 
interpreted as being due to enhanced r\ meson production 
in proton-neutron scattering [f|. For the higher energies 
accessible at CERN SPS and BNL RH1C (with ^/sWn 
ranging from 20 GeV to 200 GeV per nucleon in cm., 
in CERES, HELIOS/3, NA60 [J and PHENIX Q ex- 
periments) the abnormal dilepton yield has not been yet 
explained satisfactorily by known processes in hadronic 
physics [111]. 

Following Q we divide the range of dilepton invariant 
masses into high ( M > 3.2 GeV), low ( M < 1.2 GeV) 
and intermediate. The low mass region (LMR) is in turn 
divided into LMR I with M < 0.3 GeV and LMR II with 
0.3 GeV < M < 1.2 GeV. In LMR I the enhancement 
could possibly be explained by modifications of meson 
properties in nuclear medium [6j-|9[ as well as by proton- 
neutron scattering [5J. But in the LMR II the p meson, 
directly via 7T7t fusion or indirectly through Dalitz pro- 
cesses, largely dominates and the in-medium effects of a 
dropping mass and/or broadening resonance seem unable 
to explain the spectacular dilepton enhancement by a 4 
to 7 factor, depending on pr and centrality (see [Hj, 1H| 
for recent attempts). 

In this letter we propose a radically different explana- 
tion of this enhancement. We suggest that the effect may 
be a manifestation of local parity breaking (LPB) in col- 
liding nuclei due to the generation of a pseudoscalar, isos- 
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inglet or neutral isotriplet, condensate whose magnitude 
depends on the dynamics of the collision. It has been sug- 
gested that such a background could be due to the topo- 
logical charge fluctuations leading to the so-called Chiral 
Magnetic Effect (CME)Q studied by lattice QCD sim- 
ulations [l3| and seemingly detected in the STAR and 
PHENIX experiments at RHIC 14], although the issue is 
far from being settled. It might be also related to pseu- 
doscalar domain walls [l5j. However the fact that the 
observed dilepton excess is almost absent for peripheral 
collisions (where the CME should be more visible) and 
maximized in cental collisions makes us believe that it 
may be due to the ephemeral formation of a bona-fide 
thermodynamic phase where parity is broken, a possibil- 
ity that has been argued for in [161 ] . 

It has been shown in [17] that a pseudoscalar field 
slowly evolving in time drastically changes the electro- 
magnetic properties of the vacuum. In particular an en- 
ergetic photon propagating in this background may decay 
on shell into dileptons. This same mechanism extended 
to vector mesons is proposed here as the source for the 
abnormal dilepton yield in the LMR, i.e. in the range 
M < 1.2 GeV, for centrality -j- 20% and for moderate 
Pt < 1 GeV [2| ■ In this letter we will concentrate in the 
LMR II, and more specifically in the region around the 
p and U! resonant contribution. 

We shall assume that a time dependent but approx- 
imately spatially homogeneous background of a pseu- 
doscalar field a(t) is induced at the densities reached in 
heavy ion collisions and we will define a 4- vector related, 
£p ~ dfj,a, for later use. a(t) could be either isosinglet 
or isotriplet or even a mixture of the two, but detailed 
calculations will be presented for the case of isosinglet 
background only. 

The appropriate framework to describe electromag- 
netic interactions of hadrons at low energies is the Vector 
Dominance Model (VDM)@, [Lsj containing the lightest 
vector mesons po and w in the SU(2) flavor sector. We do 
not include <f> meson, as its typical mean free path ~ 40 
fm makes it insensitive to medium effects. Quark-meson 
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interactions are described by 

1 r 3 

Cint = qi^V^q; Vf, = -eA^Q + -g u u,J + 9 P Pp^, 



(1) 

il, g u ~ g p = g ~ 6. These values are 
extracted from vector meson decays. The Maxwell and 
mass terms are 



(V^a) = (Aft, Uft, p M = (po) M ) : 
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where m 2 / 



n 2 — 2q 2 f 2 ~ m, 2 , . This matrix reflects 
the VMD relations at the quark level [7|, [18| . Finally, in a 
pseudoscalar time-dependent background the Lagrangian 
contains a parity-odd Chern-Simons (CS) term 



Ccs(k) 



= ^tr C ejM Vj dkVi = ^ C ejfei ^f,o <9fcVJ, 6 , (3) 

which additionally mixes photons and vector mesons due 
to LPB. For isosinglet pseudoscalar background e 2 £ = 
| CI, and the mixing matrix reads 



N ab 




det (AT) = 0. (4) 



Remarkably, TV ~ m 2 . Simple order-of- magnitude con- 
siderations indicate that £ ~ ar _1 ~ 1 MeV, taking the 
time of formation of pseudoscalar condensate r = 1 fm 
and the value of condensate of order of f n . 

For isotriplet pseudoscalar background e% = 3£r 3 , 
and the corresponding CS matrix takes the form 
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det (AT) = 0. (5) 



The VMD coefficients in (jl]),© are obtained from the 
anomalous Wess-Zumino action |19| and related to the 
phenomenology of radiative decays of vector mesons [20| . 
The ratios of matrix elements for isotriplet condensate in 
([5]) are in direct agreement with the experimental decay 
constants for the processes ttq — > 77, oj — > 7To7, po — > 



19] and for the decay lo — > iririr [2l| taken from 



7Tn 

22]. Likewise the elements in ((4]) can be, in principle, 
estimated from the decays rj — > 77, 77' — > 77, u — » 
VI i Pq ~ * VI- However, phenomenologically there ex- 
ists a strong 773 — rjo mixing effect finally resolved in the 



SU(3) flavor scheme [23]. Only the ratio of the decay 
widths uj — > 777, po — » 777 is little sensitive to the mix- 
ing and confirms the off-diagonal elements of (|4]) . In this 
letter we ignore the above mixing and restrict ourselves 
toSU(2). 

The mass-shell equations for vector mesons read 



K^ U V 1 — 
^ab y v,b — u 



k v V vJb = 0, 

,2t — 2n 



K"" ee g» v (k 2 \ - m 2 ) - km - ie^ pa (pKN, (6) 

selecting out three physical polarizations for massive vec- 
tor fields. In fact, these three polarizations contribute 
into the vector field propagators as they couple to con- 
served fermion currents. The longitudinal polarization 
is orthogonal to k^ and to the CS vector £ M 
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for k 2 > 0. The mass of this state remains undistorted. 
The transversal (circular) polarizations e± on the other 
hand satisfy 



(k 2 l - m 2 ± V(C ' k) 2 - ( 2 k 2 N) 



J* 
e±. 



(8) 



As previously mentioned we restrict ourselves to an 
isosinglet pseudoscalar background a(t). The spectrum 
can be found after the simultaneous diagonalization of 
matrices m 2 , AT and particularizing to the case C/j — 

(C,o, 0,0) 
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namely 
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Thus in the case of isosinglet pseudoscalar background 
the massless photons are not distorted when mixed with 
massive vector mesons. In turn, massive vector mesons 
split into three polarizations with masses m v _ < m v L < 
rriy + . This splitting unambiguously signifies local par- 
ity breaking as well as breaking of Lorentz invariance 
due to the time-dependent background. For large enough 
|fc| > 10e 2 m v /9g 2 ( ~ my/360^ vector meson states with 
negative polarization become tachyons. However their 
group velocity remains less than the light velocity 24] 
provided that C < 20e 2 my/9g 2 ~ my/180 « 4.3 MeV. 
For higher values of £ the vacuum state becomes unsta- 
ble, namely, polarization effects give an imaginary part 
for the vacuum energy. Note that the position of reso- 
nance poles for ± polarized mesons is moving with wave 
vector |fc| and therefore they reveal themselves as "giant" 



resonances. The enlargement of the resonant region po- 
tentially leads to a substantial enhancement of their con- 
tribution to dilepton production away from their nominal 
vacuum resonance position. 

The production rate of dileptons pairs mediated by p 
mesons takes a form similar to the one given in [7j but 
with modified propagators due to LPB, according to our 
previous discussion 
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Finite lepton mass corrections are not shown in (11 1 1) but 
included in the fits. For u) mesons a similar expression 
is used but without the two pion threshold, characteris- 
tic of the dominant coupling of the p to hadronic matter 
via two-pion fusion. A simple thermal average has been 
included but we remark that the temperature T is an ef- 
fective one that may in fact depend on the range of M, 
Pt and centrality. The coefficients c pi c u parametrize in 
an effective way the total cross-sections for vector me- 
son creation. Because they are not known with precision 
in the present setting, particularly their off-shell values, 
the relative weights are used as free parameters in the 
hadronic 'cocktail' [H, HE HH- This is the procedure used 
both by NA60 and PHENIX and we shall follow it here 
too. The usual 'cocktail' contains weights normalized to 
the peripheral collisions result (roughly agreeing with ex- 
isting pp and p- nucleus data). For semi-central and cen- 
tral collisions, particularly at low pt the p/co ratio needs 
to be enhanced by a factor 1.6 in the case of NA60 Q or 
approximately 1.8 in PHENIX [27|. 

A simple thermal average is appropriate only for cen- 
tral collisions and moderate values of pr- No serious 
attempt will be made here to extrapolate to peripheral 
processes, but an appropriate mixture of distorted and 
non-distorded vector mesons should be adequate. Like- 
wise to include a proper description of the pt spectrum 
a more realistic statistical description would be needed. 
Finally, for a proper comparison with the experimen- 
tal results our expressions would have to be propagated 
through the experimental acceptances. Nevertheless the 
general features of the LPB-induced modifications should 
be already visible in a simplified description. 

NA60 has obtained accurate results for the p spectral 
function[l| by measuring the p + p~ spectrum with un- 
precedent precision and by carefully subtracting the con- 
tributions from the 'cocktail' except the p itself. The 
corresponding data for central collisions is shown in Fig. 
1 along with the contribution from the 'cocktail' scaled by 
the 1.6 factor previously indicated and using m p = 750 
MeV. The 'cocktail' p can be obtained from (fTTj) by set- 
ting £ = and adjusting the effective temperature to 
T ~ 300 MeV [26]. However, even after optimization 
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FIG. 1: The p meson contribution into the dilepton produc- 
tion is shown for parity symmetric nuclear matter £ = and 
for local parity breaking with £ = 2 MeV compared to the 
NA60 measurement for central collisons and all pr- Both 
contributions are normalized to produce the same on-shell 
cross-section for p meson production, which is fitted to the 
experiment. The effective temperature is 300 MeV. 



of the constant to reproduce central p production, the 
agreement of the data with the 'cocktail' p is poor, reflect- 
ing the long standing problem of the insufficient dilepton 
yield. In the same figure we show, after adjusting the 
scale to reproduce the p peak, the result for C — 2 MeV. 
Clearly the agreement is much better, particularly on 
the right of the p peak. On the left of the resonace there 
is a noticeable discrepancy for M < 650 MeV approxi- 
mately; experimental points seem to be shifted upwards 
by a small constant amount with respect to the predic- 
tion from LPB. It is revealing that the upper kinematical 
limit for the Dalitz process u — > p + p~ir° is 643 MeV. 

Important contributions to dilepton production at 
lower invariant masses are the Dalitz processes ttq — > 
7e + e~,?7 — > 7? + ^, nearly saturating the produc- 
tion in the M < 300 MeV region, and cj -> 1+1^tt°. 
The latter has a partial lepton width nearly identical to 
the one of p —> and thus expected to show a sim- 
ilar behavior. The Kroll-Wada formula[28| includes the 
contribution of vector mesons to the previous Dalitz pro- 
cesses and it remains valid in the case of LPB provided 
that we replace the vector meson masses by the values 
in (fTTTj) according to the intermediate meson polarization 
(L, ±). We have checked that this contribution shows an 
enhancement but this along with all other hadronic pro- 
cesses relevant for dilepton production will be discussed 
in a separate publication. Indeed NA60 has been able to 
itentify the individual contributions from r\ and uj Dalitz 
decays and found an enhancement, particularly for the 
latter po|. 

Note that the p mean free path (1.3 fm) is much 
shorter than the expected size of the hadronic gas fireball 
(L = 5 10 fm). The life-time of the fireball is compa- 



rable tfb — L. The time spent by a ultra- relativistic 
resonance of width T inside the fireball is related to its 
mean free path ~ 1/r in such a medium. Therefore 
when Ttfb C 1 a tangible suppression of the resonant 
enhancement 1/r 2 arises. Around the distorted reso- 
nance peak a crude estimation of this suppression gives 
~ r FB T 2 in relative units. This suppression affects all 
long lived particles such as the ones entering the Dalitz 
processes for the LMR I, and it is also relevant to uj me- 
son decays away from the vacuum peak. It is also the 
ultimate reason why it is unnecessary to include the <f> 
vector resonance in the discussion. In our results, unless 
stated otherwise, we have not included any modification 
of the ui propagator due to LPB. 

Moving now to the PHENIX e+e - data@, [H| we ob- 
serve that it has a much poorer precision and a numerical 
subtraction of the 'cocktail' to determine different indi- 
vidual contributions is meaningless. For low values of pt 
and for central collisions the effective temperature quoted 
by PHENIX lies in the region 100 ~ 150 MeV for the 
range of invariant masses under consideration. Based on 
the standard 'cocktail' we apply the normalization factor 
1.8 already mentioned in order to approximately account 
for the correct p/uj ratio and use Eq. (fTT|) with £ ^ 
to include the effects of LPB. Fig. 2 shows the predicted 
e + e~ yield for z — 1 which is close to the optimal fit. 
The precision of the data does not allow for such clear 
cut conclusions as in the case of NA60 but it is clear that 
LPB noticeably improves the agreement to the data as 
compared to the 'cocktail'. 

We would like to emphasize the extraordinary simplic- 
ity of the approach presented here. The fits presented use 
the values (effective temperatures, normalizations, etc. ) 
quoted by the experiment themselves. The only free pa- 
rameter is £, which is expected to depend on the char- 
acteristics of the collision. It should also be said clearly 
that the presence of LPB does not preclude other many 
body or in-medium corrections Q, as long as they do not 
represent double counting. 

We now summarize the signatures and outline possible 
searches of local parity breaking. 

Polarization: Dileptons produced for values of the in- 
variant mass above and below the p + uj pole are pre- 
dominantly of opposite circular polarizations. Thus one 
could search for a asymmetries among longitudinal and 
transverse polarization for different M in event-by-event 
measurements. These measurements may reveal in an 
unambiguous way the existence of parity violation . 

Distorted photons: At low energies massive vectors 
dominate the amplitudes but photons should also show a 
distortion induced by LPB. Those with '+' polarization 
exhibit different momentum thresholds ~ 4m 2 /C to show 
a resonant behavior for different dilepton species (l7j. 
Note that finite size suppresion is relevant for photons. 

The nature of the condensate: Mixing of photons with 
vector mesons is sensitive to isospin of pseudoscalar con- 
densate and therefore the fraction of distorted photon 
decays helps to disentangle its isospin content. 
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FIG. 2: The meson contribution to dilepton production is 
shown for parity symmetric nuclear matter £ = (discontinu- 
ous line) and for local parity breaking with £ = 1, 2 MeV (solid 
line and dotted line, respectively) compared to the PHENIX 
measurements for minimum bias, all pr events. The ( = line 
is just the 'cocktail' contribution as quoted by PHENIX [29^. 
The C 7^ results from enhancing the p/uj ratio by a factor 
1.8 without changing the oj normalization w.r.t. the hadronic 
'cocktail'. The uj contribution itself includes a 20% of LPB 
contribution = 1) and 80% pure 'cocktail' = 0) to par- 
tially account for medium effects. A reasonable fit is obtained 
for T — 110 MeV (compatible with PHENIX estimates for low 
p T ) and ( = 1 MeV. Data for M < 650 MeV is not shown 
as the contribution from the Dalitz processes has not been 
included. Inset: the rho spectral function for f = 1 MeV is 
compared to the cocktail one with the same assumptions. 

To summarize, in a time-dependent pseudoscalar back- 
ground massless photons of '+' polarization and massive 
vector mesons behave as giant resonances after averaging 
over thermal distribution. For an isosinglet pseudoscalar 
background in the framework of VMD only the massive 
vector mesons p and U) propagators are distorted due to 
mixing. We have computed their contribution and found 
that they naturally tend to produce an overabundance of 
dilepton pairs in the p + uj resonance region. The modi- 
fied p spectral function shows features very similar to the 
ones measured by NA60 in dimuon events. The PHENIX 
data is much better described by this mechanism than by 
adjusting the standard hadronic 'cocktail' in the p + uj 
resonance region. At lower invariant masses the Dalitz 
processes ir° — > r ye + e~,rj — > -fl + l~ and w — > l + l~ir° sat- 
urate the hadronic contribution to e + e~ production and 
are enhanced by the LPB-induced modifications on the 
vector meson propagators; detailed results arc deferred 
to a future publication. The only free parameter is £, 
characterizing the time variation of the pseudoscalar con- 
densate. A good fit to the data is obtained for natural 
values of The possibility that the LPB condensate is 
an isotriplet or an admixture of isotriplet and isosinglet 
has been discussed. Experimental signals of the manifes- 
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tation of LPB in heavy ion collisions have been suggested. 

Thus local parity breaking seems capable of explaining 
in a natural way the PHENIX/CERES/NA60 'anomaly' 
and searching for its manifestation in dilepton production 
represents an interesting challenge for some of the LHC 
collaborations. 

We acknowledge the financial support from projects 



FPA2007-66665, 2009SGR502, CPAN (Consolider 
CSD2007-00042) and FLAVIANET. A.& V. Andrianov 
are supported also by Grants RFBR 09-02-00073-a and 
10-02-00881-a and by SPbSU grant 11.0.64.2010. We 
thank D. d'Enterria and J. Casalderrey for discussions 
and useful remarks. 



[1] M. Masera, [HELIOS/3 Collaboration], Nucl. Phys. A 
590, 103c (1995); G. Agakichiev et al. [CERES Collabo- 
ration], Eur. Phys. J. C4, 231 (1998); Phys. Lett. B 690, 
(2010); R. Arnaldi et al.[NA60 Collaboration], Phys. 
Rev. Lett. 96, 162302 (2006). 

[2] A. Adare et al.[PHENIX Collaboration], Phys. Rev. C81, 
034911 (2010). 

[3] K. O. Lapidus and V. M. Emel'yanov, Phys. Part. Nucl. 

40, 29 (2009). 
[4] I.Tserruya, 0903.0415 [nucl-ex]. 

[5] K. O. Lapidus [HADES Collaboration], Phys. Atom. 

Nucl. 73, 985 (2010). 
[6] G. E. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 

(1991). 

[7] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25, 1 (2000); 

W. Liu, R. Rapp, Nucl. Phys. A796, 101 (2007); H. van 

Hees and R. Rapp, Nucl. Phys. A806, 339 (2008). 
[8] W. Cassing and E. Bratkovskaya, Nucl. Phys. A807, 214 

(2008); E. L. Bratkovskaya, W. Cassing and O. Linnyk, 

Phys. Lett. B670, 428 (2009). 
[9] K. Dusling and I. Zahed, Nucl. Phys. A825, 212 (2009). 
[10] H. van Hees and R. Rapp, Phys. Rev. Lett. 97, 102301 

(2006). 

[11] J. Ruppert et al, Phys. Rev. Lett. 100, 162301 (2008). 

[12] D. Kharzeev, R. D. Pisarski and M. H. G. Tytgat, Phys. 
Rev. Lett. 81, 512 (1998); K. Buckley, T. Fugleberg, 
and A. Zhitnitsky, Phys. Rev. Lett. 84, 4814 (2000); D. 
Kharzeev, Phys.Lett. B633, 260 (2006); D. E. Kharzeev, 
L. D. McLerran and H. J. Warringa, Nucl. Phys. A803, 
227 (2008). 

[13] P. Buividovich, M. Chernodub, E. Luschevskaya and M. 

Polikarpov, Phys. Rev. D 80, 054503 (2009). 
[14] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 

103, 251601 (2009); S. A. Voloshin, J. Phys. Conf. Ser. 



230, 012021 (2010). 
[15] A. Gorsky and M. B. Voloshin, 1006.5423 [hep-th]. 
[16] A. A. Andrianov and D. Espriu, Phys. Lett. B 663, 450 

(2008); A. A. Andrianov, V. A. Andrianov and D. Espriu, 

Phys. Lett. B 678, 416 (2009). 
[17] A. A. Andrianov, D. Espriu, P. Giacconi and R. Soldati, 

JHEP 0909, 057 (2009); A. A. Andrianov, D. Espriu, F. 

Mescia and A. Renau, Phys. Lett. B 684, 101 (2010) . 
[18] J.J. Sakurai, Ann. of Physics 11, 1 (1960); Currents 

and Mesons, (The University of Chicago Press, Chicago 

1969). 

[19] N. Kaiser and U.-G. Meissner, Nucl. Phys. A 519, 671 
(1990); E. Truhlik, J. Smejkal and F.C. Khanna, Nuclear 
Physics A 689, 741 (2001). 

[20] O. Dumbrais et al., Nucl. Phys. B 216, 277 (1983). 

[21] F. Klingl, N. Kaiser and W. Weise, Z. Phys. A356, 193 
(1996). 

[22] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 
075021 (2010). 

[23] T. Feldmann, P. Kroll and B. Stech, Phys. Rev. D58, 

114006 (1998) ; Y. N. Klopot, A. G. Oganesian and O. 

V. Tervaev. larXiv:0911. 0180 [hep-ph]. 
[24] J. Alfaro, A. A. Andrianov, M. Cambiaso, P. Giacconi 

and R. Soldati, Int. J. Mod. Phys. A 25, 3271 (2010). 
[25] R. Arnaldi et al [NA60 Collaboration], Phys. Lett. B 677, 

260 (2009) 

[26] R. Arnaldi et al. [NA60 Collaboration], Eur. Phys. J. C 

61, 711 (2009). 
[27] V.L. Rykov, private communication. 
[28] L. G. Landsberg, Phys. Rep. 128, 301 (1985). 
[29] Numerical data extracted from 

http://www.phenix.bnl.gov/phenix/WWW/info/data/ 

ppg088_data.html/fig26.txt 



